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Abstract. It is well known that Clifford (geometric) algebra offers a geo- 
metric interpretation for square roots of —1 in the form of blades that 
square to minus 1. This extends to a geometric interpretation of quater- 
nions as the side face bivectors of a unit cube. Systematic research has 
been done [S^ on the biquaternion roots of —1, abandoning the restric- 
tion to blades. Biquaternions are isomorphic to the Clifford (geometric) 
algebra C£{3,0) of R^. Further research on general algebras C£{p,q) has 
explicitly derived the geometric roots of —1 for p + q < 4 ,17L The cur- 
rent research abandons this dimension limit and uses the Clifford algebra 
to matrix algebra isomorphisms in order to algebraically characterize the 
continuous manifolds of square roots of —1 found in the different types 
of Clifford algebras, depending on the type of associated ring (R, H, R^, 
H^, or C). At the end of the paper explicit computer generated tables of 
representative square roots of —1 are given for all Clifford algebras with 
n = 5, 7, and s = 3 (mod 4) with the associated ring C. This includes, 
e.g., C£{Q,5) important in Clifford analysis, and C£(4, 1) which in ap- 
plications is at the foundation of conformal geometric algebra. All these 
roots of —1 are immediately useful in the construction of new types of 
geometric Clifford Fourier transformations. 
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1. Introduction 

The young London Goldsmid professor of applied mathematics W. K. Clifford 
created his geometric algebraai^ in 1878 inspired by the works of Hamilton on 

^In his original publication (Sj Clifford first used the term geometric algebras. Subsequently 
in mathematics the new term Clifford algebras |24l has become the proper mathematical 
term. For emphasizing the geometric nature of the algebra, some researchers continue [S] 
I13II14| to use the original term geometric algebra(s). 
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quaternions and by Grassmann's exterior algebra. Grassmann invented the 
antisymmetric outer product of vectors, that regards the oriented parallelo- 
gram area spanned by two vectors as a new type of number, commonly called 
bivector. The bivector represents its own plane, because outer products with 
vectors in the plane vanish. In three dimensions the outer product of three 
linearly independent vectors defines a so-called trivector with the magnitude 
of the volume of the parallelepiped spanned by the vectors. Its orientation 
(sign) depends on the handedness of the three vectors. 

In the Clifi'ord algebra [H] of the three bivector side faces of a 
unit cube {6162,6263,6361} oriented along the three coordinate directions 
{61, 62, 63} correspond to the three quaternion units i, j, and k. Like quater- 
nions, these three bivectors square to minus one and generate the rotations 
in their respective planes. 

Beyond that Clifford algebra allows to extend complex numbers to 
higher dimensions [4lfT4j and systematically generalize our knowledge of com- 
plex numbers, holomorphic functions and quaternions into the realm of Clif- 
ford analysis. It has found rich applications in symbolic computation, physics, 
robotics, computer graphics, etc. [3111131111111]. Since bivectors and trivec- 
tors in the Clifford algebras of Euclidean vector spaces square to minus one, 
we can use them to create new geometric kernels for Fourier transformations. 
This leads to a large variety of new Fourier transformations, which all deserve 
to be studied in their own right [ilinillSllTlIinilinilllllSHinilSI]. 

In our current research we will treat square roots of —1 in Clifford alge- 
bras C£{p, q) of both Euclidean (positive definite metric) and non-Euclidean 
(indefinite metric) non-degenerate vector spaces, K" ~ R"'" and M.P''^, re- 
spectively. We know from Einstein's special theory of relativity that non- 
Euclidean vector spaces are of fundamental importance in nature |12) . They 
are further, e.g., used in computer vision and robotics [5] and for general 
algebraic solutions to contact problems [33]. Therefore this chapter is about 
characterizing square roots of —1 in all Clifford algebras C£{p,q), extending 
previous limited research on C£{3, 0) in [32] and C£{p, q),n^ p+q < 4 in JJ, . 
The manifolds of square roots of —1 in C£{p, q), n = p+q = 2, compare Table 
1 of [T7] , are visualized in Fig. [1] 

First, we introduce necessary background knowledge of Clifford algebras 
and matrix ring isomorphisms and explain in more detail how we will char- 
acterize and classify the square roots of —1 in Clifford algebras in Section [H 
Next, we treat section by section (in Sections [3] to [T]) the square roots of — 1 
in Clifford algebras which are isomorphic to matrix algebras with associated 
rings M, H, R^, H^, and C, respectively. The term associated means that the 
isomorphic matrices will only have matrix elements from the associated ring. 
The square roots of —1 in Section [7] with associated ring C are of particular 
interest, because of the existence of classes of exceptional square roots of — 1, 
which all include a nontrivial term in the central element of the respective 
algebra different from the identity. Section [7] therefore includes a detailed 
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Figure 1. Manifolds of square roots / of -1 in C£{2,0) 
(left), C£{1, 1) (center), and C£{0, 2) = H (right). The square 
roots are f = a + hiei + 6262 + /3ei2, with a, 61, 62, /3 G K, 
a = 0, and /S^ = 6f e| + 6|ef + e\el. 



discussion of all classes of square roots of —1 in the algebras C£(4, 1), the iso- 
morphic C^(0,5), and in C£(7,0). Finally, we add appendix \K\ with tables of 
square roots of —1 for all Clifford algebras with n — 5, 7, and s = 3 (mod 4). 
The square roots of —1 in Section [7] and in Appendix 1X1 were all computed 
with the Maple package CLIFFORD ^ , as explained in Appendix [BJ 



2. Background and problem formulation 

Let C£{p,q) be the algebra (associative with unit 1) generated over R by 
p + q elements (with k — 1,2, . . . ,p + q) with the relations = 1 ii k < p, 
ej. = —\ \i k > p and ChCk + CkCh = whenever h ^ k, see We set the 
vector space dimension n ~ p + q and the signature s = p — q. This algebra 
has dimension 2", and its even subalgebra C£o{p,q) has dimension 2"'"^ (if 
n > 0). We are concerned with square roots of —1 contained in C£{p,q) or 
C£o{p,q). If the dimension of C£{p,q) or, C£o{p,q) is < 2, it is isomorphic 
to R ^ C£{0, 0), R2 ^ C£{1, 0),OT C£{0, 1), and it is clear that there is 
no square root of —1 in R and R^ = R x R, and that there are two squares 
roots i and — i in C. Therefore we only consider algebras of dimension > 4. 
Square roots of —1 have been computed explicitly in [351 for C£{3,0), and in 
[It] for algebras of dimensions 2" < 16. 

An algebra C£{p, q) or C£o{p, q) of dimension > 4 is isomorphic to one 
of the five matrix algebras: M{2d,R), M{d,B.), M{2d,M.'^), Mid,^^) or 
Ai{2d,C). The integer d depends on n. According to the parity of n, it is 
either 2("-2)/2 q^. 2(«-3)/2 fo^. Ci{p,q), and, either 2("-4)/2 or 2("-3)/2 f^j. 
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C^o (p, '?)■ The associated ring (either M, H, M^, H^, or C) depends on s in 
this wajo: 



s mod 8 
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Therefore we shaU answer this question: What can we say about the square 
roots of —1 in an algebra A that is isomorphic to A^(2(i, R), H), 
A4(2d,R2), M{d,¥?), or, A4(2d,C)? They constitute an algebraic subman- 
ifold in A] how many connected component^ (for the usual topology) does 
it contain? Which are their dimensions? This submanifold is invariant by 
the action of the group Inn(^) of inner automorphism^ of A, i.e. for ev- 
ery r e ^, = -1 ^ /(r)^ = -1 V/ G Inn(^). The orbits of Inn(^) are 
called conjugacy classeso; how many conjugacy classes are there in this sub- 
manifold? If the associated ring is R^ or or C, the group Aut(^) of all 
automorphisms of A is larger than Inn(^), and the action of Aut(^) in this 
submanifold shall also be described. 

We recall some properties of A that do not depend on the associated 
ring. The group Inn(yi) contains as many connected components as the group 
G(^) of invertible elements in A. We recall that this assertion is true for 
A^(2c?, R) but not for J\4{2d + 1,R) which is not one of the relevant ma- 
trix algebras. If / is an element of A, let Cent(/) be the centralizer of /, 
that is, the subalgebra of all 5 G ^ such that fa = gf- The conjugacy 
class of / contains as many connected component^ as G(^) if (and only if) 
Cent(/) Pi G(^) is contained in the neutralO connected component of G(^), 



^Compare chapter 16 on matrix representations and periodicity of 8, as well as Table 1 on 
p. 217 of [24l . 

•^Two points are in the same connected component of a manifold, if they can be joined by 
a continuous path inside the manifold under consideration. (This applies to all topological 
spaces satisfying the property that each neighborhood of any point contains a neighborhood 
in which every pair of points can always be joined by a continuous path.) 
■^An inner automorphism / of ^ is defined as / : ^ ^ ji,f{x) = a~^xa,\/x £ ji, with 
given fixed a £ The composition of two inner automorphisms g{f{x)) = xab = 

(afe)~^x(afe) is again an inner automorphism. With this operation the inner automorphisms 
form the group Inn(yt), compare \35'. 

^The conjugacy class (similarity class) of a given r G A,r^ = —1 is {/(r) ; / € Inn(^)}, 
compare |34J. Conjugation is transitive, because the composition of inner automorphisms 
is again an inner automorphism. 

^According to the general theory of groups acting on sets, the conjugacy class (as a topo- 
logical space) of a square root / of —1 is isomorphic to the quotient of G(»4) and Cent(/) 
(the subgroup of stability of /). Quotient means here the set of left handed classes modulo 
the subgroup. If the subgroup is contained in the neutral connected component of G(,4), 
then the number of connected components is the same in the quotient as in G(v4). See also 
0. 

"^Neutral means to be connected to the identity element of A.. 
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and the dimension of its conjugacy class is 

dim(yl) -dim(Cent(/)). 



(2.1) 



Note that for invertible g e Cent(/) we have g~^fg — f ■ 

Besides, let Z(^) be the center of A, and let [A, A] be the subspace 
spanned by all [f,g] = Ls — gf ■ In all cases A is the direct sum of 7j{A) 
and [A, A]. For example^ Z(7W(2d, K)) = {al | a £ E} and Z(A^(2d,C)) = 
{cl I c G C}. If the associated ring is K or H (that is for even n), then 
'L[A) is canonically isomorphic to M, and from the projection A — >■ Z(^) we 
derive a linear form Seal : — > M. When the associated rin^ is or 
or C, then Z(^) is spanned by 1 (the unit matri^l^) and some element to 
such that up' = ±1. Thus, we get two linear forms Seal and Spec such that 
Scal(/)1 + Spec{f )uj is the projection of / in Z{A) for every f £ A. Instead 
of u) we may use — w and replace Spec with —Spec. The following assertion 
holds for every f € A: The trace of each multiplicatiort"! g^fgovg^gf 
is equal to the product 



The word "trace" (when nothing more is specified) means a matrix trace 
in R, which is the sum of its diagonal elements. For example, the matrix 
M e 7W(2d,M) with elements m^; e M, 1 < fc, Z < 2d has the trace tr(M) = 



We shall prove that in all cases Scal(/) = for every square root of — 1 
in A. Then, we may distinguish ordinary square roots of —1, and exceptional 
ones. In all cases the ordinary square roots of —1 constitute a uniqueflj conju- 
gacy class of dimension dim(^)/2 which has as many connected components 
as G[A)^ and they satisfy the equality Spec(/) = if the associated ring is 
or or C. The exceptional square roots of —1 only exisi[ilif A = M{2d, C). 
In A4(2d, C) there are 2c? conjugacy classes of exceptional square roots of —1, 
each one characterized by an equality Spec(/) = k/d with ±fc e {1, 2, . . . , d} 
[see Section [7], and their dimensions are < dim(^)/2 [see eqn. (|7.5p ]. For 
instance, w (mentioned above) and — w are central square roots of —1 in 
A^(2d, C) which constitute two conjugacy classes of dimension 0. Obviously, 
Spec(w) = 1. 



*A matrix algebra based proof is e.g., given in I33| . 

^This is the case for n (and s) odd. Then the pseudoscalar ui G Ci{p, q) is also in Z(C£(p, q)). 
^''The number 1 denotes the unit of the Clifford algebra A, whereas the bold face 1 denotes 
the unit of the isomorphic matrix algebra M . 
^^These multiplications are bilinear over the center of yl. 

^■^Let A be an algebra A4(m, K) where K is a division ring. Thus two elements / and g of 
A induce K-linear endomorphisms /' and g' on K™; if K is not commutative, K operates 
on K'" on the right side. The matrices / and g are conjugate (or similar) if and only if 
there are two K-bases Bi and B2 of K."^ such that /' operates on Bi in the same way as 
g' operates on B2. This theorem allows us to recognize that in all cases but the last one 
(with exceptional square roots of —1), two square roots of —1 are always conjugate. 
^•^The pseudoscalars of Clifford algebras whose isomorphic matrix algebra has ring or 
square to = +1. 



tr{f )=dim{A) Scal(/). 



(2.2) 



Efe=l"^fefe [21]. 



6 



E. Hitzer, J. Hclmstetter and R. Ablamowicz 



For symbolic computer algebra systems (CAS), like MAPLE, there exist 
Clifford algebra packages, e.g., CLIFFORD [3], which can compute idempo- 
tents [2J and square roots of —1. This will be of especial interest for the 
exceptional square roots of —1 in M{2d, C). 

Regarding a square root r of — 1, a Clifford algebra is the direct sum 
of the subspaces Cent(r) (all elements that commute with r) and the skew- 
centralizer SCent(r) (all elements that anticommute with r). Every Clifford 
algebra multivector has a unique split by this Lemma. 

Lemma 2.1. Every multivector A € C£{p, q) has, with respect to a square root 
r G Ci{p,q) of —1, i.e., r~^ = — r, the unique decomposition 

A± ^^{A±r-^Ar), A = A+ + A^, A+r = rA+, A^r = -rA^. (2.3) 

Proof. For A G C£{p, q) and a square root r G C£{p, q) of —1, we compute 

A±r^ ^{A±r^^Ar)r ^ i(Ar ± r"M(-l)) ^{rr^^Ar ±rA) 

= ±r^{A±r-'^Ar). 



For example, in Clifford algebras C£{n, 0) [2^ of dimensions n = 2 mod 4, 
Cent(r) is the even subalgebra Cloin, 0) for the unit pseudoscalar r, and the 
subspace C£i{n, 0) spanned by all fc-vectors of odd degree k, is SCent(r). The 
most interesting case is A^(2c?, C), where a whole range of conjugacy classes 
becomes available. These results will therefore be particularly relevant for 
constructing Clifford Fourier transformations using the square roots of —1. 



3. Square roots of -1 in M{2d, R) 

Here A = M{2d,R), whence dim{A) = {2d)^ = Ad^ . The group G{A) has 
two connected components determined by the inequalities det{g) > and 
det(g) < 0. 

For the case d = 1 we have, e.g., the algebra C£{2,0) isomorphic to 
M{2,R). The basis {1,61,62,612} of C£(2,0) is mapped to 

{(0 o)'(o o^)}' 

The general element a + bici + 6262 + (3ei2 G C£{2, 0) is thus mapped to 

in A^(2,R). Every element f of A = A^(2d, M) is treated as an ]R- linear 
endoniorphism of = M^'^. Thus, its scalar component and its trace (|2.2|) 
arc related as follows: tr(/) — 2c?Scal(/). If / is a square root of —1, it turns 
V into a vector space over C (if the complex number i operates like fonV). 
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If (ei, 62, . . . , ed) is a C-basis of V, then (ei, /(ei), 62, /(e2), . . . , e^, /(ed)) is 
a R-basis of and the 2d x 2(i matrix of / in this basis is 




d 



Consequently all square roots of — 1 in are conjugate. The centralizer of a 
square root / of —1 is the algebra of all C-linear endomorphisms g oiV (since 
i operates like / on V). Therefore, the C-dimension of Cent(/) is cP and its 
M-dimension is 2cP. Finally, the dimension (|2.1I) of the conjugacy class of / is 
dim{A) - dim(Cent(/)) = 4^^ - 2d'^ = 2d'^ = dim(y^)/2. The two connected 
components of G(,4) are determined by the sign of the determinant. Because 
of the next lemma, the M-determinant of every element of Cent(/) is > 
0. Therefore, the intersection Cent(/) P| G(^) is contained in the neutral 
connected component of G(^) and, consequently, the conjugacy class of / 
has two connected components like G(^). Because of the next lemma, the K- 
trace of / vanishes (indeed its C-trace is di, because / is the multiplication 
by the scalar i: f{v) = iv for all v) whence Scal(/) = 0. This equality is 
corroborated by the matrix written above. 

We conclude that the square roots of —1 constitute one conjugacy class 
with two connected components of dimension dim(^)/2 contained in the 
hyperplane defined by the equation 

Scal(/) = 0. (3.3) 

Before stating the lemma that here is so helpful, we show what happens 
in the easiest case d = 1. The square roots of —1 in A^(2,R) are the real 
matrices 

il -»)(* -»)-(«' + '>")i-^^ M 

hence a? + be = —1, a relation between a, 6, c which is equivalent to {h — 
c)^ = (& + c)^ + 4a^ + 4 {h - > 4 b ~ c > 2 (one component) 
or c — > 2 (second component). Thus, we recognize the two connected 
components of square roots of —1: The inequality b > c + 2 holds in one 
connected component, and the inequality c > & + 2 in the other one, compare 
Fig.H 

In terms of C€(2, 0) coefScients with b-c = -(-/3+&i) = 2/3, 
we get the two component conditions simply as 

/? > 1 (one component), f3 < —1 (second component). (3-5) 

Rotations {det{g) = 1) leave the pseudoscalar /3ei2 invariant (and thus pre- 
serve the two connected components of square roots of —1), but reflections 
{det{g') = —1) change its sign /3ei2 — > — /3ei2 (thus interchanging the two 
components) . 

Because of the previous argument involving a complex structure on 
the real space V, we conversely consider the complex space with its 
structure of vector space over R. If (ei, 62, . . . , e^) is a C-basis of C'', then 
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Figure 2. Two components of square roots of —1 in (2,1 



(ei, iei, 62, i&ii • ■ • , e^, ierf) is a M-basis. Let g be a C-linear endomorphism of 
(i.e., a complex d x d matrix), let trc(5) and detc{g) be the trace and 
determinant of g in C, and trR((/) and detR((7) its trace and determinant for 
the real structure of C^. 

Example. For d = 1 an endomorphism of is given by a complex number 
g — a + ib, a, b £ R. Its matrix representation is according to 

with ( : "1 = (a^ - 62) 



-b 
a 



1 
1 



2ab 



-1 

1 



Then we have trc(g) = a+ib, trg 



2a 



(3.6) 

23?(trc(g)) and detc(.9) 



a + ib, deti! 



= a2 + 62 = |detc(ff)l2 >0. 



Lemma 3.1. For every C-linear endomorphism g we can write tr-g^g) = 
23?(trc(5)) and detm{g) = | detc{g)\^ > 0. 

Proof. There is a C-basis in which the C-matrix of g is triangular [then 
detc(fl') is the product of the entries of g on the main diagonal]. We get the 
R-matrix of g in the derived R-basis by replacing every entry a + of the 

C-matrix with the elementary matrix . The conclusion soon follows. 

The fact that the determinant of a block triangular matrix is the product of 
the determinants of the blocks on the main diagonal is used. □ 



4. Square roots of -1 in M{2d,M.'^) 

Here A = M{2d,M.^) = M{2d,R) x X(2d,R), whence dim(^) = Sd^. 
The group G(^) has foui0 connected components. Every element (/, /') G 
A (with /, /' e Ai{2d,R)) has a determinant in which is obviously 

^*In general, the number of connected components of G(.4) is two if A = A1(m,K), and 
one if A = A4(m, C) or ^ = A^(m,H), because in all cases every matrix can be joined 
by a continuous path to a diagonal matrix with entries 1 or —1. When an algebra A is 
a direct product of two algebras B and C, then G(^) is the direct product of G{B) and 
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(det(/), det(/')), and the four connected components of G(^) are determined 
by the signs of the two components of detR2(/, /'). 

The lowest dimensional example {d = 1) is C£{2,1) isomorphic to 
A4{2,M.'^). Here the pseudoscalar w = 6123 has square oj^ = +1. The cen- 
ter of the algebra is and includes the idempotents e± = (l±aj)/2, 
ej_ = e±, e+e_ = e_e+ — 0. The basis of the algebra can thus be written 
as {e+, eie+, 626+, 6126+, e_, eie_, 626-, 6126-}, where the first (and the last) 
four elements form a basis of the subalgebra C£{2, 0) isomorphic to M{2, R). 
In terms of matrices we have the identity matrix (1,1) representing the scalar 
part, the idempotent matrices (1, 0), (0, 1), and the w matrix (1, —1), with 1 
the unit matrix of M{2,R). 

The square roots of (—1,-1) in A are pairs of two square roots of — 1 
in M{2d,M.). Consequently they constitute a unique conjugacy class with 
four connected components of dimension 4d^ — dim(^)/2. This number 
can be obtained in two ways. First, since every element (/, /') G A (with 
/, /' G Ai{2d,R)) has twice the dimension of the components / G A^(2d, M) 
of Section [31 we get the component dimension 2 • 2d^ = 4c?^. Second, the cen- 
tralizer Cent(/, /') has twice the dimension of Cent(/) of A4(2d, M), therefore 
dim(^) - Cent(/, /') = - Ad"^ = Ad^ . In the above example for d = 1 
the four components are characterized according to p. 51) by the values of the 
coefficients of /3ei2e+ and /3'ei2e- as 



Cl : 


/3> 1, 




C2 : 


/3> 1, 


/?' < -1, 


ca : 


/3<-l, 


/?'>!, 


C4 : 


/3<-l, 


P' < -1- 



For every (/,/') e A we can with ^ write tr(/) + tr(/') = 2dScal(/,/') 
and 

tr(/)-tr(/') = 2dSpec(/,/') if c. = (l,-l); (4.2) 

whence Scal(/, /') ~ Spec(/, /') = if {f,f') is a square root of (—1, — 1), 
compare p.3p . 

The group Aut(^) is larger than Inn(^), because it contains the swap 
automorphism (/,/') 1— > (/',/) which maps the central element w to — w, 
and interchanges the two idempotents e+ and e_. The group Aut{A) has 
eight connected components which permute the four connected components 
of the submanifold of square roots of (—1,-1). The permutations induced 
by Inn(^) are the permutations of the Klein group. For example for d = 1 of 



G(C), and the number of connected components of G(.A) is the product of the numbers of 
connected components of G{B) and G(C). 
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(|4.ip we get the following Inn(Al(2, R^)) permutations 

det(g) > 0, det(5') > : identity, 

det(.g) > 0, det(5') < : (ci, C2), (cg, C4), 

det(g) < 0, det(5') > : (ci, C3), (c2, C4), 

det(.g)<0, det(5')<0: (ci, C4), (c2, C3). (4.3) 

Beside the identity permutation, Inn(^) gives the three permutations that 
permute two elements and also the other two ones. 
The automorphisms outside Inn(^) are 

{fJ')^{9f'g'\9'f9''') for some {g,g')eG{A). (4.4) 

If det(sr) and det(g') have opposite signs, it is easy to realize that this auto- 
morphism induces a circular permutation on the four connected components 
of square roots of (—1, —1): If det(5) and det(g') have the same sign, this 
automorphism leaves globally invariant two connected components, and per- 
mutes the other two ones. For example, for d = 1 the automorphisms (14. 4p 
outside Inn(^) permute the components (|4.ip of square roots of (—1, —1) in 



M{2, 



■) as follows 
det(5) > 0, 
det(5) > 0, 
det(5) < 0, 
det(5) < 0, 



det(.g') > 
det(g') < 
det(g') > 
det(g') < 



(Cl), (C2,C3), (C4), 
Ci ^ C2 ^ C4 -J> C3 
Cl ^ C3 ^ C4 -!> C2 
(ci,C4), (C2), (C3). 



Cl, 
Cl, 



(4.5) 



Consequently, the quotient of the group Aut(^) by its neutral connected 
component is isomorphic to the group of isometrics of a square in a Euclidean 
plane. 



5. Square roots of —1 in Ai{d,M) 

Let us first consider the easiest case d = 1, when ^ = H, e.g., of C£(0, 2). The 
square roots of —1 in H are the quaternions ai + bj + cij with a^-f = 1. 

They constitute a compact and connected manifold of dimension 2. Every 
square root / of —1 is conjugate with i, i.e., there exists v G M : v~^fv = 
i ^ fv = vi. If we set v — ^fi + 1 — a + bij — cj + 1 we have 

fv = -ft + f = f + i = ifi-i) + l)i = vi. 

V is invertible, except when / = —i. But i is conjugate with —i because 
ij = hence, by transitivity / is also conjugate with — i. 

Here A = M{d,M), whence dim(A) = Ad'^. The ring H is the algebra 
over R generated by two elements i and j such that = = — 1 and 
ji = —ij. We identify C with the subalgebra generated b\P^ i alone. 

The group G(^) has only one connected component. We shall soon 
prove that every square root of —1 in A is conjugate with i\. Therefore, the 



This choice is usual and convenient. 
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submanifold of square roots of —1 is a conjugacy class, and it is connected. 
The centralizer of il in ^ is the subalgebra of all matrices with entries in C. 
The C-dimension of Cent(zl) is cP , its R-dimension is 2(P, and, consequently, 
the dimension ()2.1|) of the submanifold of square roots of —1 is AcP — 2(P = 
2d^ = dim(A)/2. 

Here V — is treated as a (unitary) module over H on the right side: 
The product of a line vector 'w = {xi,X2, ■ ■ ■ ,Xd) € V hy y G M is *v y = 
{xiy, X2y, • ■ ■ , Xdy)- Thus, every f G A determines an H-linear endomorphism 
of V: The matrix / multiplies the column vector v = *{xi,X2, ■ • ■ , Xd) on the 
left side v ^ fv. Since C is a subring of H, V is also a vector space of dimen- 
sion 2d over C. The scalar i always operates on the right side (like every scalar 
in H). If (ei, 62, ... , e^) is an H-basis of V , then (ei, eij, 62, 62^, . . . , e^, e^j) 
is a C-basis of V . Let / be a square root of —1, then the eigenvalues of / in 
C are +z or —i. If we treat V as a. 2d vector space over C, it is the direct 
(C-linear) sum of the eigenspaces 

V+ ^{veV\ f{v) = vi} and V- ^ {veV \ f{v) = -vi}, (5.1) 

representing f as a, 2d x 2d C-matrix w.r.t. the C-basis of V, with C-scalar 
eigenvalues (multiplied from the right): X± = ±i. 

Since ij = —ji, the multiplication v 1— )■ vj permutes V'^ and ^ as 
f{v) — ±vi is mapped to f{v)j = ±vij = T('fj)*- Therefore, if (ei, 62, ... , e^) 
is a C-basis of , then {eij, 623, ■ ■ ■ , erj) is a C-basis of V~ , consequently 
(ei, eij, 62, e2j, . . . , e^, e^j) is a C-basis of V , and (ei, 62, ... , er=d) is an H- 
basis of V . Since / by /(e^) = gfci for fc = 1, 2, . . . , d operates on the H-basis 
(ei, 62, ... , Ed) in the same way as il on the natural H-basis of V , we conclude 
that / and «1 are conjugate. 

Besides, Scal(il) = because 2il = G [A,A\, thus il ^ Z{A). 

WhenceEl 

Scal(/) — for every square root of — 1. (5.2) 

These results are easily verified in the above example of d = 1 when 
^ = H. 

6. Square roots of -1 in M{d,M^) 

Here, A = M{d,W^) = M{d,m) x M(d,H), whence dim(A) = . The 
group G(^) has only one connected component (see Footnote [Tl)) . 

The square roots of (—1,-1) in A are pairs of two square roots of — 1 
in A^(rf, H). Consequently, they constitute a unique conjugacy class which is 
connected and its dimension is 2 x 2d^ — Ad^ = dim(^)/2. 

For every (/, /') G A we can write Scal(/) + Scal(/') = 2 Scal(/, /') 
and, similarly to (|4.2p . 

Scal(/)-Scal(/') = 2Spec(/,/') if c^ = (l,-l); (6.1) 



^''Compare the definition of Scal(/) in Section[2l remembering that in the current section 
the associated ring is H. 
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whence Scal(/, /') — Spec(/, /') = if {f,f') is a square root of (—1,-1), 
compare with (|5.2p . 

The group Aut(^) has twcFI connected components; the neutral com- 
ponent is Inn(^), and the other component contains the swap automorphism 

(/,/') ^ (/',/)• 

The simplest example is d ~ 1, A = H^, where we have the identity 
pair (1, 1) representing the scalar part, the idempotents (1, 0), (0, 1), and w 
as the pair (1, —1). 

^ = is isomorphic to C£{0,3). The pseudoscalar u = 6123 has the 
square oj^ = +1. The center of the algebra is {l,aj}, and includes the idem- 
potents e± — i(l±a;), e^. ~ e±, e+e_ = £-£+ = 0. The basis of the algebra 
can thus be written as {e_|., eie+, 626+, 6126+, e_, eie_, 626-, ei2e_} where the 
first (and the last) four elements form a basis of the subalgebra C£{0, 2) iso- 
morphic to H. 



7. Square roots of -1 in M{2d, C) 

The lowest dimensional example for d = 1 is the Pauli matrix algebra A ~ 
A4(2,C) isomorphic to the geometric algebra C£(3,0) of the 3D Euclidean 
space and C^(l,2). The C£{3,0) vectors 61,62,63 correspond one-to-one to 
the Pauli matrices 

^i=(j J), ^2=(^ 7), a3^(^j \y (7.1) 

with critT2 = ia^ = (t ^ ■] ■ The element uj = ci 172(73 = il represents the 



central pseudoscalar 6123 of C£(3, 0) with square = —1. The Pauli algebra 
has the following idempotents 

ei^af^l, eo = i(l + CT3), £-1=0. (7.2) 
The idempotents correspond via 

/ = z(2e-l), (7.3) 

to the square roots of —1: 

»),/.^.<'3^(; °,),/_,. 

where by complex conjugation /_i = /i. Let the idempotent Eq = ^{1 — (73) 
correspond to the matrix /q = —473. We observe that /o is conjugate to 
/o = "'f /o^i = '^i'^2 = fo using = cti but /i is not conjugate to 
/_i. Therefore, only /i,/o,/-i lead to three distinct conjugacy classes of 
square roots of —1 in A^(2,C). Compare Appendix IB] for the corresponding 
computations with CLIFFORD for Maple. 



^'^ Compare Footnote 1141 
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In general, if = M{2d, C), then dim{A) = Sd^. The group G{A) has 
one connected component. The square roots of —1 in ^ are in bijection with 
the idempotents e [2] according to (|7.3p . Accordin^ll to (|7.3p and its inverse 
e — i(l — if) the square root of —1 with Spec(/_) — k/d = ~l, i.e. k ^ —d 
(see below), always corresponds to the trival idempotent e_ = 0, and the 
square root of —1 with Spec(/+) = k/d = +1, k = +d, corresponds to the 
identity idempotent e+ = 1. 

If / is a square root of —1, then V ~ C^"* is the direct sum of the 
eigenspaceS associated with the eigenvalues i and —i. There is an integer 
k such that the dimensions of the eigenspaces are respectively d + k and 
d — k. Moreover, —d<k<d. Two square roots of —1 are conjugate if and 
only if they give the same integer k. Then, all elements of Cent(/) consist of 
diagonal block matrices with 2 square blocks of {d + k) x (d + k) matrices 
and (d — k) X {d — k) matrices. Therefore, the C-dimension of Cent(/) is 
{d + fc)^ + (d— fc)^. Hence the M-dimension (|2.ip of the conjugacy class of /: 

- 2{d + kf - 2{d - kf = 4(^2 _ fc2), (7,5) 

Also, from the equality tr(/) = (d + k)i — (d — k)i = 2ki we deduce that 
Scal(/) = and that Spec(/) = {2ki)/{2di) = k/d if w = il (whence 
tr(a;) = 2di). 

As announced on page[5l we consider that a square root of —1 is ordinary 
if the associated integer k vanishes, and that it is exceptional if fc ^ . Thus 
the following assertion is true in all cases: the ordinary square roots of —1 in 
A constitute one conjugacy class of dimension dim(^)/2 which has as many 
connected components as G(„4), and the equality Spec(/) — holds for every 
ordinary square root of —1 when the linear form Spec exists. All conjugacy 
classes of exceptional square roots of —1 have a dimension < dim(^)/2. 

All square roots of —1 in A^(2d, C) constitute (2d+l) conjugacy classej^ 
which are also the connected components of the submanifold of square roots 
of —1 because of the equality Spec(/) = k/d, which is conjugacy class spe- 
cific. 

When A = A4(2d, C), the group Aut(^) is larger than Inn(^) since 
it contains the complex conjugation (that maps every entry of a matrix to 
the conjugate complex number). It is clear that the class of ordinary square 
roots of —1 is invariant by complex conjugation. But the class associated 



^°On the other hand it is clear that complex conjugation always leads to /_ = /+, where 
the overbar means complex conjugation in M{2d,C) and Clifford conjugation in the iso- 
morphic Clifford algebra C£{p,q). So either the trivial idempotent e_ = is included in 
the bijection 117.31 1 of idempotents and square roots of —1, or alternatively the square root 
of —1 with Spec(/_) = — 1 is obtained from /_ = /+. 

^^The following theorem is sufficient for a matrix / in Ai{m,K.), if K is a (commutative) 
field. The matrix / is diagonalizable if and only if P{f) = for some polynomial P that 
has only simple roots, all of them in the field K. (This implies that P is a multiple of 
the minimal polynomial, but we do not need to know whether P is or is not the minimal 
polynomial). 

■^''Two conjugate (similar) matrices have the same eigenvalues and the same trace. This 
suffices to recognize that 2d+ I conjugacy classes are obtained. 
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with an integer k other than is mapped by complex conjugation to the 
class associated with —k. In particular the complex conjugation maps the 
class {uj} (associated with fc = d) to the class {— w} associated with k — — d. 

All these observations can easily verified for the above example of d = 1 
of the Pauli matrix algebra A = Ai{2,C). For d — 2 we have the isomor- 
phism of A = A4(4,C) with C£{0,5), C£(2,3) and C^(4,l). While C£{0,5) 
is important in Clifford analysis, C£{4:, 1) is both the geometric algebra of 
the Lorentz space R*^^ and the conformal geometric algebra of 3D Euclidean 
geometry. Its set of square roots of —1 is therefore of particular practical 
interest. 



Example. Let C^(4, 1) = A where A = X(4,C) for d 
1-vectors can be represented^ by the following matrices: 



2. The C£(4,l) 



ei = 



/I 








-1 








\0 






64 






-1 



/o 



1 
Vo 






V 





-1 



62 = 



1 



vo 

-1 



0/ 





1 



1 



, 63 





-1 



/o 

i 


Vo 





—i 
0/ 



o\ 

1 



0/ 



(7.6) 



We find five conjugacy classes of roots ft of —1 in C€(4, 1) for k e {0, ±1, ±2}: 
four exceptional and one ordinary. Since fk is a root of p{t) — + 1 which 
factors over C into (t — i){t + i), the minimal polynomial mk{t) of fk is 
one of the following: t — i, t + i, or (t ~ i){t + i). Respectively, there are 
three classes of characteristic polynomial Afe(i) of the matrix in yV4(4,C) 
which corresponds to fk, namely, {t — i)"^, (t + iY, and {t — i)"''-{t + i)^^, where 
ni + n2 ~ 2d — A and ni = d + k = 2 + k, 712 — d — k = 2 — k. As predicted 
by the above discussion, the ordinary root corresponds to fc = whereas the 
exceptional roots correspond to A: 7^ 0. 

1. For k — 2, we have A2(t) — {t ~ i)**, m2{t) — t — i, and so 7^2 = 
diag(z, i, i, i) which in the above representation ()7.6p corresponds to the 
non-trivial central element /2 = w = 612345. Clearly, Spec(/2) = 1 = ^; 
Scal(/2) ~ 0; the C-dimension of the centralizer Cent(/2) is 16; and the 
R-dimension of the conjugacy class of /2 is zero as it contains only /2 
since /2 £ Z{A). Thus, the M-dimension of the class is again zero in 
agreement with (17. 5p . 

2. For k = -2, we have A_2(0 = (t + i)", m_2(t) = t + i, and J'_2 = 
diag(— i, —i, —i, —i) which corresponds to the central element /_2 = 



^^For the computations of this example in the Maple package CLIFFORD we have used the 
identification i = 623 ■ Yet the results obtained for the square roots of — 1 are independent 
of this setting (we can alternatively use, e.g., i = £12345, or the imaginary unit i £ C), as 
can easily be checked for /i of 1)7. 7p . fo of 1)7. 8p and /-i of l|7.9p by only assuming the 
standard Clifford product rules for ei to 65. 
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= -612345- Again, Spec(/_2) = -1 = ^; Scal(/__2) = 0; the C- 
dimension of the centralizer Cent(/_2) is 16 and the conjugacy class of 
f-2 contains only /_2 since /_2 6 Z{A). Thus, the R-dimension of the 
class is again zero in agreement with ()7.5p . 
3. For k ^ ±2, we consider three subcases when fc = 1, fc = 0, and k = —1. 
When fc = 1, then Ai(i) = {t - if{t + i) and mi(t) = [t - i){t + i). 
Then the root J^i = diag(i,i,i, —i) corresponds to 

/l = 2^^'^^ + '^123 - 62345 + 612345)- (7-7) 

Note that Spec(/i) = ^ = ^ so /i is an exceptional root of —1. 
When fc = 0, then Ao(i) = {t-ifit + if and TOo(i) = {t-i){t + i). Thus 
the root of —1 in this case is J-q = diag(i, i, — i, — i) which corresponds 
to just 

/o = 6123- (7.8) 

Note that Spec(/o) = thus /o = 6123 is an ordinary root of —1. 
When fc = -1, then A_i(i) = {t-i){t + if and r7i_i(t) = {t-i){t + i). 
Then, the root of —1 in this case is T-i = diag(z, — i, — i, — i) which 
corresponds to 

/-I = 2 ("^23 + ei23 + 62345 - 612345)- (7-9) 

Since Scal(/_i) = — ^ = ^, we gather that /_i is an exceptional root. 

As expected, we can also see that the roots lo and — w are re- 
lated via the grade involution whereas /i = — /-i where ~ denotes the 
reversion in C£(4, 1). 

Example. Let C^(0, 5) ^ A where A = >[(4,C) for d = 2. The C£(0,5) 
1-vectors can be representee^ by the following matrices: 



/o 


-1 













—i 













/ 


—i 





\ 


1 















—i 




































-1 


, 62 













—i 




63 = 












i 







1 













J 








v 








-i/ 








/o 





-1 








/ 


— i 

























1 












I 








(7-10) 




64 




1 











, 65 




-i 








1 












-1 





0^ 













0^ 









Like for C£(4, 1), we have five conjugacy classes of the roots fk of —1 in 
C£(0,5) for fc G {0, ±1,±2}: four exceptional and one ordinary. Using the 



^^For the computations of this example in the Maple package CLIFFORD we have used the 
identification i = 63. Yet the results obtained for the square roots of —1 are independent 
of this setting (we can alternatively use, e.g., i = £12345, or the imaginary unit i £ C), as 
can easily be checked for /i of H7.11| l. /o of l|7.12p and /_i of 1)7. 13p by only assuming the 
standard Clifford product rules for ei to 65. 
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same notation as in Example [71 we find the following representatives of the 
conjugacy classes. 

1. For k — 2,we have A2{t) — (t—i)'^, m2{t) = t—i, and J-2 = diag(i, i, i, i) 
which in the above representation (I7.10p corresponds to the non-trivial 
central element /2 = = £12345. Then, Spec(/2) = 1 = |; Scal(/2) — 0; 
the C-dimension of the centralizer Cent(/2) is 16; and the R-dimension 
of the conjugacy class of /2 is zero as it contains only /2 since /2 G Z{A). 
Thus, the M-dimension of the class is again zero in agreement with ()7.5p . 

2. For k = -2, we have A„2(0 = {t + i)'^, m_2(i) = t + i, and J'_2 = 
diag(— — i, —i, —i) which corresponds to the central element /-2= — 
w = -612345- Again, Spec(/_2) = -1 = f ; Scal(/_2) = 0; the C- 
dimension of the centralizer Cent(/_2) is 16 and the conjugacy class of 
f-2 contains only /_2 since /_2 G Z{A). Thus, the R-dimension of the 
class is again zero in agreement with ()7.5p . 

3. For k ^ ±2, we consider three subcases when A; = 1, fc = 0, and k = —1. 
When fc = 1, then Ai(t) = (t - if{t + i) and mi{t) = {t - i){t + i). 
Then the root J-i = diag(«,i,«, — «) corresponds to 

/i ^(es + ei2 + 645 + 612345)- (7-11) 

Since Spec(/i) = | = /i is an exceptional root of —1. 
When k = Q, then Ao(i) = {t-if{t + if and mo(t) = {t-i){t + i). Thus 
the root of —1 is this case is J-q ~ diag(i, i, — i, —i) which corresponds 
to just 

/o - 645- (7.12) 

Note that Spec(/o) = thus /o = 645 is an ordinary root of —1. 
When k = -1, then A„i(t) = {t-i){t + if and r7i_i(<) = {t-i){t + i). 
Then, the root of —1 in this case is T-i = diag(z, — z, — i, — z) which 
corresponds to 

/-I = ^(^63 + 612 + 645 - 612345)- (7-13) 

Since Scal(/_i) = — ^ = ^, we gather that /_i is an exceptional root. 

Again we can see that the roots /2 and f-2 are related via the 
grade involution whereas /i = — /-i where ~ denotes the reversion in 
C^(0,5). 

Example. Let C^(7, 0) = A where A = A4(8,C) for d = 4. We have nine 
conjugacy classes of roots fk of —1 for k G {0, ±1, ±2 ± 3 ± 4}. Since is 
a root of a polynomial p{t) ~t^ + l which factors over C into [t — i){t + i), 
its minimal polynomial m(t) will be one of the following: t — i, t + i, or 
{t-i){t + i) =t^ + l. 

Respectively, each conjugacy class is characterized by a characteristic 
polynomial Ak{t) of the matrix Mk G A^(8, C) which represents fk - Namely, 
we have 

Afe(^) = (^-^)"H^ + ^)"^ 
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where ni + n2 = 2d = 8 and ni = d + k = 4: + k and n2 = d— k = A — k. 
The ordinary root of —1 corresponds to fc = whereas the exceptional roots 
correspond to fc 7^ 0. 

1. When fc = 4, we have A4(i) — [t — i)^, m4{t) = t — i, and J-4 = 

8 

diag(i, ■ ■ ■ ,i) which in the representation used by CLIFFORD W cor- 
responds to the non-trivial central element f4=uj = 61234567- Clearly, 
Spec(/4) = 1 = ^; Scal(/4) = 0; the C-dimension of the centralizer 
Cent(/4) is 64; and the K-dimension of the conjugacy class of /4 is zero 
since /4 S Z{A). Thus, the M-dimension of the class is again zero in 
agreement with ()7.5p . 

2. When k = -4, we have A_4(t) = (t + i)^, m_4(t) = t + i, and 



= diag(— i, . . . , —i) which corresponds to /_4 = — u; = —61234567- 
Again, Spec(/_4) ~ —1 = 2', Scal(/_4) — 0; the C-dimension of the 
centralizer Cent(/) is 64 and the conjugacy class of /_4 contains only 
/_4 since /_4 £ Z{A). Thus, the M-dimension of the class is again zero 
in agreement with (j7.5l) . 
3. When k 7^ ±4, we consider seven subcases when k = ±3, k — ±2, 
k = ±1, and fc = 0. 

When fc = 3, then A^it) = (t - + i) and m3(i) = (t - i){t + i). 

7 

Then the root J-3 = diag(i, . . . ,i, —i) corresponds to 

/s = ^(^'23 - 645 + 667 - 6123 + 6145 - 6167 + 6234567 + 36i234567)- (7-14) 

Since Spec(/3) = | = ^, /s is an exceptional root of —1. 

When fc 2, then A2{t) = (t ~ if{t + if and m2(t) = (i - -f i). 

6 

Then the root J-2 = diag(i, . . . , i, — i, — i) corresponds to 

h = ^(667 - 645 - 6123 + 61234567)- (7-15) 

Since Spec(/2) = 5 = /2 is also an exceptional root. 

When fc = 1, then Ai{t) = (t - if{t + if and mi{t) = (i - i){t + i). 

5 

Then the root Ti = diag(i, . . . , — i, — i, — i) corresponds to 

/l = ^(^23 - 645 + 3667 - 6123 + 6145 + 6167 " 6234567 + 6i234567)- (7-16) 

Since Spec(/i) = |; = |, /i is another exceptional root. 

When fc = 0, then Ao(t) = (t - iY{t + if and moit) = (i - i){t + i). 

Then the root J-q = diag(i, i, i, i, —i, —i, —i, —i) corresponds to 

/o = ^(623 - 645 + 667 - 6234567)- (7-17) 

Since Spec(/o) = = ^, we see that /o is an ordinary root of —1. 
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When k = -1, then A_i(i) {t-if{t + iy^ and TO_i(t) = {t~i){t + i). 

5 

Then the root J-^i = diag(i, i, i, —i, . . . , —i) corresponds to 

/-I = ^(^23 — 645 + 3e67 + 6123 — 6145 — 6157 — 6234567 " 61234567)- (7-18) 

Thus, Spec(/_i) = — i = 2 and so /_i is another exceptional root. 
When k = -2, then A_2(0 = + and m_2(i) = + 

6 

Then the root J- ^2 = diag(i, — i, . . . , — i) corresponds to 

f-2 = ^(667 - 645 + 6123 " 61234567)- (7-19) 

Since Spec(/_2) = —\ = see that /_2 is also an exceptional root. 

When k = -3, then A_3(i) = {t-i){t + iy and m_3(t) = {t-i){t + i). 

7 

Then the root J-'„3 = diag(i, — i, . . . , — i) corresponds to 

/-3 = ^(^23 - 645 + 667 + 6123 - 6145 + 6167 + 6234567 " 3ei234567)- (7.20) 

Again, Spec(/_3) = ~| = ^ and so /_3 is another exceptional root 
of -1. 

As expected, we can also see that the roots lo and —uj are related 
via the reversion whereas = —f-3, /2 = ^/-2, fi = ^.f-i where 
denotes the conjugation in C£{7,0). 



8. Conclusions 

We proved that in all cases Scal(/) = for every square root of —1 in ^ 
isomorphic to Ci{p^q). We distinguished ordinary square roots of —1, and 
exceptional ones. 

In all cases the ordinary square roots / of — 1 constitute a unique con- 
jugacy class of dimension dim(^)/2 which has as many connected compo- 
nents as the group G(,4) of invertible elements in A. Furthermore, we have 
Spec(/) — (zero pseudoscalar part) if the associated ring is M^, H^, or C. 
The exceptional square roots of —1 only exist if ^ = A4{2d,C) (see Sec- 
tion [T]). 

For A = Ai{2d,M) of Section [31 the centralizer and the conjugacy class 
of a square root / of — 1 both have R-dimension 2d^ with two connected 
components, pictured in Fig. [2] for d — 1. 

For A = X(2(i,M2) = M(2d,R) x M{2d,R) of Section H the square 
roots of (—1,-1) are pairs of two square roots of —1 in A4{2d,M.). They 
constitute a unique conjugacy class with four connected components, each of 
dimension Ad'^. Regarding the four connected components, the group Inn(^) 
induces the permutations of the Klein group whereas the quotient group 
Aut(^) /lnn{A) is isomorphic to the group of isometrics of a Euclidean square 
in 2D. 
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For A = ^A{d, H) of Section [5l the submanifold of the square roots / 
of —1 is a single connected conjugacy class of M-dimension equal to the 
M-dimcnsion of the centralizer of every /. The easiest example is HI itself for 
d=l. 

For A = M{d,m^) = 7W(2d,H) x M{2d,M.) of Section El the square 
roots of (—1,-1) are pairs of two square roots {f,f') of —1 in A^(2d,]HI) 
and constitute a unique connected conjugacy class of R-dimension id"^. The 
group Aut(^) has two connected components: the neutral component Inn(^) 
connected to the identity and the second component containing the swap 
automorphism (/, /') i— > (/', /). The simplest case for d = 1 is H'^ isomorphic 
to C£(0,3). 

For A — A4{2d, C) of Section [3 the square roots of —1 are in bijection 
to the idempotents. First, the ordinary square roots of —1 (with k = 0) con- 
stitute a conjugacy class of R-dimension of a single connected component 
which is invariant under Aut(^). Second, there are 2d conjugacy classes of 
exceptional square roots of —1, each composed of a single connected compo- 
nent, characterized by equality Spec(/) = k/d (the pseudoscalar coefficient) 
with ±fc G {1, 2, . . . , d}, and their R-dimensions are 4((i^ — fc^). The group 
Aut(^) includes conjugation of the pseudoscalar w h- > —u which maps the 
conjugacy class associated with k to the class associated with —k. The sim- 
plest case for d = 1 is the Pauli matrix algebra isomorphic to the geometric 
algebra C£(3, 0) of 3D Euclidean space R^, and to complex biquaternions 

m- 

Section [7] includes explicit examples for d — 2: C^(4, 1) and C£(0,5), 
and for d — A: C£{7,0). Appendix 1X1 summarizes the square roots of —1 
in all C£{p,q) = A4{2d,C) for d = 1,2,4. Appendix IbI contains details on 
how square roots of —1 can be computed using the package CLIFFORD for 
Maple. 

Among the many possible applications of this research, the possibility 
of new integral transformations in Clifford analysis is very promising. This 
field thus obtains essential algebraic information, which can e.g., be used to 
create steerable transformations, which may be steerable within a connected 
component of a submanifold of square roots of —1. 

Appendix A. Summary of roots of —1 in Ci{p, q) = Ai{2d, C) 
for d = 1,2,4 

In this appendix we summarize roots of —1 for Clifford algebras C£{p,q) = 
M(2d, C) for d = 1, 2, 4. These roots have been computed with CLIFFORD [3]. 
Maple [30] worksheets written to derive these roots are posted at [18]. 
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k 


fk 




1 


OJ = ei23 


{t-if 





623 


{t-i){t + i) 


-1 




(1 + rf 



Table 1. Square roots of -1 in Ci{3,0) ^ M{2,C), d=l 



k 


fk 


A,Xt) 


2 


W = 612345 




1 


^(623 + 6i23 — 62345 + 612345) 


{t-if{t + i) 





6123 


it-i)\t + if 


-1 


^(623 + 6123 + 62345 — 612345) 


{t-i){t + if 


-2 


— W = —612345 





Table 2. Square roots of -1 in C^(4, 1) ^ A^(4,C), d = 2 



k 


fk 


Afc(i) 


2 


W = 612345 




1 


5(63 + 612 + 645 + 612345) 


{t-i)^{t + i) 





645 


{t-i)^{t + i)^ 


-1 


|(-63 + 612 + 645 - 612345) 


{t-i){t + i)^ 


-2 


— W = —612345 


{t+iy 



Table 3. Square roots of -1 in C^(0,5) ^ A4(4,C), d = 2 



k 


fk 


Afc(i) 


2 


^ — 612345 


(t-z)4 


1 


5(63 + 6134 + 6235 +w) 


(<-z)3(t + .) 





6134 


(t-z)2(t + z)2 


-1 


5(-(;! + (U:!l + '^2:i--, --■) 


(i-i)(i + ,:f 


-2 


— W = —612345 


{t+if 



Table 4. Square roots of -1 in C^(2,3) ^ A^(4,C), = 2 
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k 


/fe 


Afc(i) 


4 


^ 61234567 




3 


-6167 + 6234567 + 3cj) 


{t - iy{t + i) 


2 


1(667 - £45 - ei23 + ^) 


(t-lf{t + lf 


1 


4- f69'^ — 6d.f^ -|- 36fi7 — ei9q 4- e^A^ 

+ 6167 — ^234567 + ^) 


{t - iy{t + i)^ 





1(^23 — ^45 + 667 — 6234567) 




-1 


■T f69T — 6/if; 4- 36r7' -(- 619"? — 6i/i^ 

— 6167 — 6234567 ^ 


it - i)^(t + i)^ 


-2 


5 (667 - 645 + 6123 " 


(t- 1^(1 + if 


-3 


I: (623 - 645 + 667 + 6123 - 6145 

+ 6167 + 6234567 — 3uj) 


{t-i){t+iy 


-4 


— —61234567 




Table 5. Square roots of -1 in C£{7,0) ^ 7W(8,C), d = 4 


k 


fk 


Afc(i) 


4 


W = 61234567 




3 


3: (64 - 623 - 656 + 61237 + 6147 

+ 61567 — 623456 + 3w) 




2 


^(-623 - 656 + ei47 + w) 


{t-if{t + if 


1 


|;(-64 - 623 - 3656 " 61237 + 6l47 

+ 61567 — 623456 + 


{t-if{t + if 





5(64 + 623 + 656 + 623456) 


{t-i)\t + iy 


-1 


J (-64 - 623 - 3656 + 61237 - 6147 

— 61567 — 623456 " 


(t~if{t + if 


-2 


5 (-623 - 656 - 6147 - w) 




-3 


I; (64 — 623 — 656 — 61237 — 6147 

^ 61567 ~ 623456 — 3w) 


{t-i){t + iy 


-4 


—<-^ = —61234567 





Table 6. Square roots of -1 in C£(l, 6) ^ A^(8, C), d = 4 



Appendix B. A sample Maple worksheet 

In this appendix we show a computation of roots of —1 in C^(3, 0) in CLIF- 
FORD. Although these computations certainly can be performed by hand, 
as shown in Section [71 they illustrate how CLIFFORD can be used instead 
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k 


fk 


Afc(i) 


4 


^ — 61234567 




3 


— ^13457 — 623467 + 3a;) 


{t - iy{t + i) 


2 


^(ei45 — ^12456 — ^13457 + ^) 


{t-lf{t + if 


1 


— 613457 — 623467 + ^) 







|(64 + 612456 + 613457 + 623467) 


{t-l)\t + lf 


-1 


~7 ( — -1 — " 1 -1 1^ O /I Q /I 7 — " 3 1 OA'^R 

4 V '-^4 '-'•±40 '-^Z4D 1 '-'o4 f '-"-■iZ40D 

— 613457 — 623467 — 


it - i)^(t + i)^ 


-2 


^( — 6145 — 612456 " 613457 — u) 


(t- 1^(1 + if 


-3 


3: (64 — 6145 — 6246 — 6347 — 612456 

— 613457 — 623467 — 3a;) 


{t-i){t+iy 


-4 


— —61234567 




Table 7. Square roots of -1 in C£{3,A) = 7W(8,C), d = 4 


k 


h 




4 


OJ = ei234567 




3 


3: ( — 623 + ei23 + 62346 + 62357 — 612346 

— 612357 + 6234567 + Suj) 


{t~iy{t + i) 


2 


^(6123 — 612346 — 612357 + w) 


{t-ifa + if 


1 


3: ( — 623 + 6123 — 62346 + 62357 — 3ei2346 
— 612357 — 6234567 + 


{t-if{t + if 





§(623 + 612346 + 612357 + 6234567) 




-1 


3: ( — 623 — 6123 + 62346 — 62357 — 3612346 
— 612357 — 6234567 — 




-2 


^( — 6123 — 612346 — 612357 — Oj) 


it-iy(t + if 


-3 


■J ( — 623 — 6123 — 62346 — 62357 — 612346 

— 612357 + 6234567 — 3l^) 


(t-i){t + iy 


-4 


— = —61234567 


{t + ^y 



Table 8. Square roots of -1 in C£(5, 2) = M{8, C), d = 4 



especially when extending these computations to higher dimensionsHl To see 



^•^In showing Maple display we have edited Maple output to save space. Package asvd is 
a supplementary package written by the third author and built into CLIFFORD. The pri- 
mary purpose of asvd is to compute Singular Value Decomposition in Clifford algebras [l] . 
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the actual Maple worksheets where these computations have been performed, 
see [IB]. 

> restart :with(Clif ford) :with(linalg) :with(asvd) : 

> p,q:=3,0; ##« — selecting signature 

> B:=diag(l$p,-l$q) : ##« — defining diagonal bilinear form 

> evaKmakealiases (p+q) ) : ##<< — defining aliases 

> clibas : =cbasis (p+q) ; ##assigning basis for CI (3,0) 



p, g := 3, 

chbas ■- [Id, el, e2, e3, el2, el3, e23, el23] 



> data: =clidata() ; ##« — displaying information about Cl(3,0) 

data ■— [complex, 2, simple, + [Id, e2, e3, e23], [Id, e23], [Id, e2]] 

> MM : =matKrepr ( ) ; ##<< — displaying default matrices to generators 

Cliplus has been loaded. Definitions for type/climon and 
type/clipolynom now include &C and &C[K]. Type ?cliprod for help. 



1 





, e2 = 





1 


, e3 = 





-~e23 





-1 




1 







e23 






MM ■- [el 



Pauli algebra representation displayed in (7.1): 

> Sigma [1] :=evalm(rhs(MM[l] )) : 

> Sigma [2] : =evalm (rhs (MM [2] ) ) : 

> Sigma [3] : =evalm (rhs (MM [3] ) ) : 






1 




1 





1 





e23 



-e23 






1 


7 


-1 



We show how we represent the imaginary unit i in the field C and the diagonal 
matrix diag(i, i) : 

> ii:=e23; ##« — complex imaginary unit 

> II : =diag(ii , ii) ; ##« — diagonal matrix diag(i,i) 

II := e23 



II 



e23 
e23 



We compute matrices mi,m2, ■ ■ ■ representing each basis element in 
C^(3,0) isomorphic with C(2). Note that in our representation element 623 
in Cl{Z, 0) is used to represent the imaginary unit i. 
> for i from 1 to nops (clibas) do 

lprint('The basis element clibas [i] is represented by the following 

matrix : ' ) ; 

M[i] : =subs (Id=l .matKrepr (clibas [i] ) ) od; 
'The basis element'. Id, 'is represented by the following matrix:' 

1 



Ml := 



1 
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"The basis element', el, 'is represented by the following matrix:' 

M2 := 



1 
-1 



"The basis element', e2, 'is represented by the following matrix:' 



Ms 



1 

1 



"The basis element', e3, 'is represented by the following matrix:' 

M4 := 



-e23 
e23 



"The basis element', el2, 'is represented by the following matrix:' 

Ms := 



1 
-1 



"The basis element', el3, 'is represented by the following matrix:' 



Me 



-e23 
-e23 



'The basis element', e23, 'is represented by the following matrix:' 

M7 := 



e23 
-e23 



"The basis element', el23, 'is represented by the following matrix:' 



Ms 



e23 
e23 



We will use the procedure phi from the asvd package which gives an 
isomorphism from C(2) to C£(3, 0). This way we can find the image in C£(3, 0) 
of any complex 2x2 complex matrix A. Knowing the image of each matrix 
mi, m2, - . . , TOs in terms of the Clifford polynomials in C^(3, 0), we can easily 
find the image of A in our default spinor representation of C£(3, 0) which is 
built into CLIFFORD. 

Procedure Centralizer computes a centralizer of / with respect to the 
Clifford basis L: 
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> Centralizer:=proc(f ,L) local c,LL,m,vars,i,eq,sol; 
m:=acld(c[i] *L[i] ,i=l. .nops(L)) ; 

vars: = [seq(c[i] ,i=l. .nops(L))] ; 
eq:=clicollect(cmul(f ,m)-cmul(m,f )) ; 
if eq=0 then return L end if : 
sol:=op(clisolve(eq,vars)) ; 
m:=subs(sol,m) ; 
m:=collect (m,vars) ; 

return sort( [coeff s(m,vars)] ,bygrade) ; 
end proc : 

Procedures Seal and Spec compute the scalar and the pseudoscalar 
parts of /. 

> Scal:=proc(f ) local p: return scalarpart (f ) ; end proc: 

> Spec:=proc(f ) local N; global p,q; 
N:=p+q: 

return coeff (vectorpart (f ,N) , opCcbasis (N,N) ) ) ; 
end proc: 

The matrix idcmpotcnts in C(2) displayed in (7.2) are as follows: 

> d:=l :Eps [1] :=sigina[l] &cm sigma[l] ; 

> Eps[0] :=evalm(l/2*(l+sigma[3])); 

> Eps[-1] :=diag(0,0) ; 



Eps^, EpsQ, Eps_i :-- 



This function ff computes matrix square root of —1 corresponding to the 
matrix idempotent eps: 

> f f :=eps->evalm(II Stem (2*eps-l)); 

ff := eps — evalm(77 &cm (2 eps — 1)) 

Wc compute matrix square roots of —1 which correspond to the idempotents 
Epsi, Epso, Eps-i, and their characteristic and minimal polynomials. Note 
that in Maple the default imaginary unit is denoted by /. 

> F [1] : =f f (Eps [1] ) ; ##<< — this square root of -1 corresponds to Eps[l] 
Delta [1] :=charpoly(subs(e23=I,evalm(F[l] )) ,t) ; 

Mu[l] :=minpoly(subs(e23=I,evalm(F[l])) ,t) ; 



1 




1 










1 


5 





7 









e23 
e23 



Ai := {t - 7)^ Ml := i - 7 



> F[0] :=ff (Eps[0] ) ; ##« — this square root of -1 corresponds to Eps [0] 
Delta [0] :=charpoly(subs(e23=I,evalm(F[0])) ,t); 
Mu[0] :=minpoly(subs(e23=I,evalm(F[0])) ,t) ; 



7^0 := 



e23 
-e23 



Ao := (t-7)(t + 7), Mo:=l + f 



> F [-1] : =f f (Eps [-1] ) ; ##« — this square root of -1 corresponds to Eps[-1] 
Delta[-1] :=charpoly(subs(e23=I,evalm(F[-l])) ,t) ; 
Mu[-1] :=minpoly(subs(e23=I,evalm(F[-l] )) ,t) ; 
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F-i := 



-e23 
-e23 



A_l:=(^ + /)^ M-i:=t + I 



Now, we can find square roots of —1 in C£{3,0) which correspond to the 
matrix square roots F^i, Fo, Fi via the isomorphism (p : C£{3,0) — >■ C(2) 
realized with the procedure phi. 

First, we let repri denote element in C£{3,0) which represents the 
diagonal {2d) x {2d) with I = i on the diagonal where i"^ = —1. This element 
will replace the imaginary unit / in the minimal polynomials. 

> reprl:=phi(<iiag(l$(2*<i)) ,M) ; 

reprI := el23 

Now, we compute the corresponding square roots /i, /o, /-i in C£{3, 0). 

> f [1] :=phi(F[l] ,M) ; ##« — element in Cl(3,0) corresponding to F[l] 
cmuKf [1] ,f [1] ) ; ##« — checking that this element is a root of -1 
Mu[l] ; ##« — recalling minpoly of matrix F[l] 

subs(e23=I,evalm(subs(t=evalm(F[l]),Mu[l]))); ##«— F[l] in Mu[l] 
mu[i] :=subs(I=reprI,Mu[i] ) ; ##« — defining minpoly of f[l] 
cmul(f [i]-reprl,ld) ; ##« — verifying that f [1] satisfies mu[l] 

/i := el23 



-Id, t-I, 






/Lti := t - el23, 



> f [0] :=phi(F[0] ,M) ; ##« — element in Cl(3,0) corresponding to F[0] 
cmuKf [0] ,f [0] ) ; ##« — checking that this element is a root of -1 
Mu[0] ; ##« — recalling minpoly of matrix F[0] 

subs (e23=I,evalm( subs (t=evalm (F [0] ) ,Mu[0] ))) ; ##«— F[0] in Mu[0] 
mu[0] :=subs(I=reprI,Mu[0] ) ; ##« — defining minpoly of f [0] 
cmuKf [0]-reprI,f [0]+reprI); ##«— f [0] satisfies mu[0] 

fo := e23 



-Id, i + r, 



Mo 



1 + r 









> f [-1] :=phi(F[-l] ,M) ; ##«— element in Cl(3,0) corresponding to F[-l] 
cmuKf [-1] ,f [-1] ) ; ##« — checking that this element is a root of -1 
Mu[-1] ; ##« — recalling minpoly of matrix F[-l] 

subs(e23=I,evalm(subs(t=evalm(F[-l]) ,Mu[-l]))) ; ##«— F[-l] in Mu[-1] 
mu[-l] :=subs(I=reprI,Mu[-l] ) ; ##« — defining minpoly of f[-l] 
cmuKf [-l]+reprl,ld) ; ##«— f [-1] satisfies mu[-l] 

/-I := -el23 



-Id, t-\-I, 






H_i -.= 1+6123, 
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Functions RdimCentralizer and RdimConjugClass of d and k compute 
the real dimension of the centralizer Cent(/) and the conjugacy class of / 
(see (7.4)). 

> RdimCentralizer: = (d,k)->2*((d+k)"2+(d-k)"2) ; ##«— from the theory 

> RdimConjugClass : = (d,k)->4*(d"2-k'2) ; ##«— from the theory 

RdimCentralizer := {d, k) ^ 2 {d + kf + 2 (d - kf 
RdimConjugClass := {d, k) —¥ 4d^ — 4:k^ 

Now, we compute the centralizers of the roots and use notation d,k,ni,n2 
displayed in Examples. 
Case k = 1 : 

> d:=l:k:=l:nl:=d+k;n2:=d-k; 

Al:=diag(I$iil,-I$n2) ; ##« — this is the first matrix root of -1 



nl := 2, n2 := 0, Al := 



I 
I 



> f [1] :=phi(Al,M); cmuKf [1] ,f [1] ) ; Scal(f[l]), Spec(f[l]); 

/i := el23, -Id, 0, 1 



> LLl:=Centralizer(f [1] .clibas) ; ##« — centralizer of f[l] 
dimCentralizer : =nops (LLl) ; ##« — real dimension of centralizer of f[l] 
RdimCentralizer (d, k) ; ##« — dimension of centralizer of f [1] from theory 

evalb(dimCentralizer=RdimCentralizer (d,k) ) ; ##« — checking equality 

LLl ■- [Id, el, e2, e3, el2, el3, e23, el23] 

dimCentralizer := 8, 8, true 

Case k = : 

> d:=l:k:=0:nl:=d+k;n2:=d-k; 

A0:=diag(I$nl,-I$n2) ; ##« — this is the second matrix root of -1 



nl 



n2 := 1, AO 



I 
-I 



> f [0] :=phi(AO,M); cmul (f [0] , f [0] ) ; Scal(f[0]), Spec(f[0]); 

fo ■- e23, -Id, 0, 



> LLO:=Centralizer(f [0] , clibas) ; ##« — centralizer of f[0] 

dimCentralizer : =nops (LLO) ; ##« — real dimension of centralizer of f [0] 
RdimCentralizer (d, k) ; ##« — dimension of centralizer of f[0] from theory 
evalb(dimCentralizer=RdimCentralizer(d,k)) ; ##« — checking equality 

LLO := [Id, el, e23, el23] 

dimCentralizer := 4, 4, true 
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Case k = —1 : 

> d:=l:k:=-l:nl:=d+k;n2:=<i-k; 

Ami :=<iiag(I$nl ,-I$n2) ; ##« — this is the third matrix root of -1 



nl := 0, n2 := 2, Ami := 



-I 
-/ 



> f [-1] :=phi(Aml,M) ; cmuKf [-1] ,f [-1] ) ; Scal(f[-1]), Spec(f[-1]); 

/_i := -el23, -Id, 0, -1 



> LLml : =Centralizer (f [-1] , clibas) ; ##« — centralizer of f [-1] 

dimCentralizer :=nops(LLml) ; ##« — real dimension of centralizer of f[-l] 
RdimCentralizer(d,k) ; ##<< — dimension of centralizer of f [-1] from theory 
evalb(dimCentralizer=RdimCentralizer (d,k) ) ; ##« — checking equality 



LLml := [Id, el, e2, e3, el2, el3, e23, el23] 
dimCentralizer := 8, 8, true 



Wc summarize roots of —1 in C£(3,0): 
'F[l] '=evalm(F[l] ) ; ##«— square root of -1 in C(2) 
Mu[l]; ##«— minpoly of matrix F[l] 
'f[l]'=f[l]; ##«— square root of -1 in Cl(3,0) 
mu[l] ; ##« — minpoly of element f[l] 



e23 
e23 
fi = el23, t - el23 



t-I 



> 'F[0] '=evalm(F[0] ) ; ##«— square root of -1 in C(2) 
Mu[0]; ##«— minpoly of matrix F[0] 
'f[0]'=f[0]; ##«— square root of -1 in CI (3,0) 
mu[0] ; ##« — minpoly of element f[0] 



e23 
-e23 



1 + r 



fo = e23, 1 + t' 



> 'F[-l] '=evalm(F[-l]) ; ##«— square root of -1 in C(2) 
Mu[-1]; ##« — minpoly of matrix F[-l] 
'f [-1] '=f [-1] ; ##«— square root of -1 in Cl(3,0) 
mu[-l] ; ##« — minpoly of element f[-l] 



-e23 






-e23 



t + I 



f-i = -el23, t+el23 
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Finaly, we verify that roots /i and /_i are related via the reversion: 
> reversionCf [1] )=f [-1] ; evalbC/,) ; 

-el23 = -el23, true 
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